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Abstract 

We compute the anomalous dimensions of field strength operators TrjF^ in A/" = 4 SYM 
from an asymptotic nested Bethe ansatz to all-loop order. Starting from the exact solution 
of the one-loop problem at arbitrary L, we derive a single effective integral equation for the 
thermodynamic L oo limit of these dimensions. We also include the recently proposed 
phase factor for the S- matrix of the planar AdS/CFT system. The terms in the effective 
equation corresponding to, respectively, the nesting and the dressing are structurally very 
similar. This hints at the physical origin of the dressing phase, which we conjecture to 
arise from the hidden presence of infinitely many auxiliary Bethe roots describing a non- 
trivial "filled" structure of the theory's BPS vacuum. We finally show that the mechanism 
for creating effective nesting/dressing kernels is quite generic by also deriving the integral 
equation for the all-loop dimension of a certain one-loop so (6) singlet state. 



1 Motivation, Conclusion and Overview 



There is much evidence that planar A/" = 4 SYM theory is integrable and that its spectral 
problem is therefore exactly solvable. It was shown by Minahan and Zarembo that the 
dilatation operator in the scalar matter sector at one loop can be mapped to the Hamilto- 
nian of an integrable so (6) spin chain and hence its eigenvalues can be found with a Bethe 
ansatz pQ. This extends to the full set of operators, leading to an integrable non-compact 
nearest-neighbor super magnet [2j . The special "solvable" properties of the A/" = 4 model 
under dilatation were already hinted at by Lipatov in f3], and extend a rather generic if 
incomplete phenomenon in more general gauge theories such as QCD, as first shown by 
Belitsky, Braun, Derkachov, Korchemsky and Manashov [3]. 

The concept of factorized scattering, one of the hallmarks of integrability, can be ex- 
tended to higher loop orders 0, and to strong coupling [6], where the gauge theory is 
expected to be more suitably described by a superstring theory in a curved AdS^ x 
background. Various analyses of this topic [7j led to a set of asymptotic all-loop Bethe 
equations for the full theory [S]. It should be stressed that quantum integrability remains 
to be proven in both gauge and string theory. In particular, on the gauge side one phe- 
nomenologically finds an integrable long-range spin chain, and the all-loop factorization 
of the multi-body magnon S-matrix into two-body processes currently has to be assumed. 
Correspondingly, symmetry fixes the magnon S-matrix only up to an overall phase factor 
[HI El E] • The latter encodes our lack of understanding of the underlying microscopic inte- 
grable structur^. However, as was argued by Janik, the dressing phase may be constrained 
by invoking crossing- invarianc^ [TT]. And indeed the string S-matrix satisfies crossing to 
the known [HI [12] orders [13]. A proposal for the complete structure of the dressing phase 
has recently been made in [13] by combining Bethe ansatz techniques for the all-order 
perturbative large spin limit of Wilson twist operators [15] with conjectures by Beisert, 
Hernandez and Lopez [TB] on the full (asymptotic) structure of the string Bethe ansatz [H] . 

An independent four-loop calculation strongly supports this guess [T7]. This calculation 
scheme has been algorithmically improved and leads to agreement between field theory and 
the dressed Bethe ansatz with a margin of error of 0.001% [L8\. Incorporating this dressing 
phase into the Bethe ansatz, an integral equation for the universal scaling function f{g) 
of A/" = 4 gauge theory in the large spin limit was derived in [13]. Some of the analytic 
properties of the solution of this equation were analyzed in [IS], and in the leading strong 
coupling limit the behavior of the scaling function as predicted by string theory [20] through 
the AdS/CFT correspondence agrees with the analysis of the equation [211 |22l EHl [23] . The 
subleading order was also successfully, albeit numerically, compared in [21] to the known 
string result [23] . 

^ As shown in [3], the full S-matrix, which is proportional to the tensor product of two copies of Shastry's 
R-matrix for the Hubbard model, satisfies the Yang-Baxter equation (see also [TU]). It should be stressed, 
however, that Yang-Baxter symmetry of the two-hody S-matrix is a necessary, but not sufficient condition 
for the integrability of a given system. An infamous counter-example is the bosonic Hubbard model. 

^ Crossing symmetry should also be a consequence of - rather than an axiom for - the proper microscopic 
formulation of the planar AdS/CFT system. 
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The dressing phase of [13] takes a surprisingly complex and seemingly opaque form. 
Arguing that it should be of a fundamental nature is, to put it mildly, unconvincing. Let 
us recall its structure. The phase shift 26{u,u') = —i log a'^{u,u') stemming from the 
dressing factor a'^{u,u') when two magnons with rapidities u,u' pass each other (see [8l fll] 
for the notation) is 

/oo /'OO / 

dte^'^'e-^ / rft'e^*'"'e-^ ikd{2 gt,2 gt') - kd{2 gt\2 gt) 
oo J —oo 

(1.1) 

where the "magic kernel" reads 

roo ±11 

kd{t,t') = 8g' dt"k,{t,2gt")^^ko{2gt",t'), (1.2) 
Jo e — i 

and the symmetric kernels Kq, Ki are expressed with the help of Bessel functions as 

,n tJi(t) Jo(tO-tVo(t) Ji(tO 

This way of writing the dressing phase has a decidedly "thermodynamic" flavor. Never- 
theless, the claim of [13] is certainly that this phase should contribute to the anomalous 
dimensions of short operator^ at and beyond four-loop order, up until "wrapping order" . A 
natural explanation of its structure could come from a hidden "non-trivial" vacuum. Recall 
that this is usually the case in relativistic field theories, or else, in "physical" ground state 
configurations of magnetic systems such as the Heisenberg or Hubbard antiferromagnets. 
What we are saying is that the dressing factor indicates that the BPS states 

TrZ^, (1.5) 

do not constitute boring "ferromagnetic" reference vacua, but are rather states more akin 
to an "antiferromagnetic" , "physical" vacuum. For reasons to be understood, BPS vacuum 
polarization at weak coupling appears for the first time at four-loop order. The difference, 
in the su(2) sector, of the AdS/CFT system to the Hubbard model [2E] is then that in the 
latter this vacuum polarization does not appear. Toy models approximatively demonstrat- 
ing a similar phenomenon in the present context have appeared previously in [271 128] . 

In the current paper we will not rigorously analyze the underlying vacuum structure 
of the BPS states (11.51) . which requires to go beyond the Bethe equations of [8J. In order 
to nevertheless prove our point, we will instead focus on another "false vacuum", but 
one which may be treated with the asymptotic spectral equations in [8]. Namely, we will 
consider operators which at one loop are built from a sole field strength component JF: 

TrJ^^. (1.6) 



^ Strictly speaking the phase was found in 14J by studying the large-spin limit of twist operators. This 
is very similar to a thermodynamic limit in that it also involves infinitely many Bethe roots ^15j . It would 
be very important to check the phase at four loops for a finite, short operator without taking any limit. 
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These are of interest as they are not embedded in a rank-one subsector, and are therefore 
described by an asymptotic nested Bethe ansatz. If L is large enough, i.e. the operator 
is "long", we may trust the ansatz to arbitrary loop order. In fact, we may pass to the 
thermodynamic limit L oo, in which the dimension of the operators (11.61) becomes 



A{g) = 2L + 8Lg' H dt ^-^^^ p{t) , (1.7) 

where p{t) is essentially the Fourier-transform of the density of momentum-carrying first- 
level Bethe roots. As is common in thermodynamic situations, the non-linear nested Bethe 
equations may be turned into a system of linear integral equations for p{t) in conjunction 
with a number of further auxiliary densities describing the distribution of the higher-level 
Bethe roots. Interestingly, the auxiliary densities may be eliminated, and one ends up with 
a single, effective linear equation for the principal density p{t). It reads 

p(t) = e-2*(l + e*)[jo(2(7t) 

-AgH dt'i-^Koi2gt,2gt') + -^K^{2gt,2gt')\pit') 

-AgH J dt' {kr,{2gt,2gt') + ka{2gt,2gt')yp{t') 

(1.8) 

Here Kd is the above kernel (11. 2p summarizing the effects of the dressing phase. Interest- 
ingly, together with a few further terms in the second line of (11.81) . the net influence of the 
nesting is exerted by the kernel Kn which reads 

poo J.II 

kn{t,t')=4g' dt"k,{t,2gt")^^ko{2gt",t'). (1.9) 
Jo e + i 

We are confident that the reader's sharp eyes will spot the structural similarity between this 
nesting kernel Kn and the dressing kernel in (11.21) . In conclusion, this result suggests 
that the AdS/CFT phase factor [H] is generated by hidden, nested levels of a yet-to-be 
constructed final Bethe ansatz. 

In the remainder of this paper we will work out the asymptotic anomalous dimensions 
of the operators (11.61) and derive the above equations. The detailed calculations are rather 
technical, and proceed as follows: In section [2] we obtain the asymptotic spectrum of these 
operators, which form a long-range spin chain, by applying the Bethe equations of [Sj. In 
section [3] we solve the one-loop problem, i.e. we find the principal and auxiliary one-loop 
Bethe roots for the appropriate Dynkin diagram. This solution is exact for all finite values 
of the length L. We also find all one-loop higher conserved charges at arbitrary L. The 
heart of the paper is section |U where we work out the thermodynamic limit of our one- loop 
solution, find (thermodynamically) all higher-loop perturbations, and reduce the resulting 
system of integral equations to the single effective one in (11.81) . Finally, in section [5], we 
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study a few further thermodynamic high-energy states satisfying special filhng conditions. 
Apart from exhibiting certain intriguing transcendentahty properties, the examples provide 
further support for our picture of the origin of nesting and dressing factors. We end with 
a short outlook. As always when working with Bethe equations, great care has to be taken 
to understand the distribution of roots. The appendices contain numerous checks and 
considerations justifying the methodology and statements in the main body of the text. 

2 Excitation Scheme for Field Strength Operators 

We are interested in the pseudo- vacuum states Tr JF-^ in (11.61) . Here pseudo- vacuum refers 
to the fact that these states may serve as a reference vacuum for the one-loop Bethe ansatz 
^J. In an oscillator realization of the system, see [22], is the tensor product of L fields 
|JF) where 

l-^) = alat|0), (2.1) 

and a{ are bosonic creation operators. After imposing the trace condition, these states are 
gauge-invariant, but, in contradistinction to the superficially similar BPS states in 
(II. 5p . not protected. In fact, their anomalous dimension is very large [22] • This type of 
operator was studied at one-loop in an integrable sector of QCD in [30] . 

There exist alternative choices for the Dynkin diagrams of a superalgebra due to the 
freedom to choose each node to be either fermionic or bosonic, with the constraint that at 
least one node has to be fermionic. Changing from one diagram to another changes the 
excitation pattern of the spin chain, and sometimes even induces a change of the corre- 
sponding vacuum state. In particular, the one-loop vacuum of the standard distinguished 
Dynkin diagram of psu(2, 2|4) is precisely given by the states Tr^"^ [2]. Therefore no Bethe 
ansatz, which takes care of the diagonalization of excitations, is necessary, and the one-loop 
anomalous dimension is just the vacuum-energy. However, it is not known how to deform 
the one-loop Bethe equations of the distinguished Dynkin diagram to higher loops. The 
same is true for most other choices, and the long-range Bethe equations of [H] only appear 
to "work" for a very specific diagram. In that diagram the vacuum is Tr^^"^, and TrjF-^ is 
a highly excited state with many excitations, whose momenta have to be diagonalized in 
order to reproduce the correct energy. 

A somewhat similar pseudo- vacuum consists of a cyclic tensor product of L fermions W, 
i.e. TyIA^ , which was studied in [HI]. However, while the latter lies in a rank-one subsector 
of the symmetry group, su(l|l), the states Tr^''" are not confined to any such subsector. 
Correspondingly, in order to diagonalize the latter one needs a nested Bethe ansatz. The 
excitation pattern of Bethe roots for the higher-loop Dynkin diagram of [8] reads 

(i^i, K2, Ks, K,, i^5, Ke, Kj) = (0, 0, 2L - 3, 2L - 2, L - 1, L - 2, L - 3) , (2.2) 

where the K,y is the excitation number of the u-th node of the Dynkin diagram. Clearly 
all but the first two nodes are highly excited. 
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Figure 1: Oscillator realization and fundamental magnons. There are four adjoint com- 
plex scalars A", A?, 3^, 3^, four light-cone covariant derivatives, P, P, and eight adjoint 
fermions lA, V,U, V,U, V,U, V. The nodes of the Dynkin diagram are labeled from 1 to 7, 
starting at the south-west end and going clockwise until the north-east end is reached. 



It is interesting to understand the state Tr J-"^ in the picture of the elementary excita- 
tions of the long-range spin chain: the 8+8 magnons of the Af = 4 model. We tabulated 
in Fig. 1 all "fundamental magnons" in a way which allows to read of the excitation pat- 
tern in the higher-loop Dynkin diagram with ease. We also included the magnon-creation 
operators in the oscillator picture, cf. [29j. The su(2|2) ©su(2|2) invariant S-matrix acts, 
respectively, on the rows and columns of Fig. 1. Exciting the central 4-th node times 
inserts X bosons into the BPS vacuum Z^. The K4 Bethe roots of the central node 
parametrize the momenta of these magnons. Next, exciting times the fifth node converts 

of the X bosons into U fermions. The process continues until all magnons contained 
in a given state are created. It follows from (12. 2p that the fields JF in Tr JF^ are not ele- 
mentary. Instead, each is essentially equivalent to a composite of two fermions U and 
V minus one field Z. Let us explain this in more detail for a single JF. The two fermions 
and the background field Z are expressed in terms of the oscillators as 

U = 4^110) "P = al4|0) Z = 44|0). (2.3) 
and form a JF, cf (12. ip . while inserting another Z: 

UV= 4A,A4,2\^)i®\^)2 ^ = al^ial_i 424^2 |0)i® 10)2. (2.4) 

We see that each JF corresponds to two magnon excitations lA V. Removing L background 
fields Z does not change the excitation number, but reduces the length (as well as the 
R-charge) by L units. Correspondingly, the state TrjF^ has, at one loop, 2L excitations 
living on a lattice of L sites. Beyond one loop, this length is not conserved. 
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In order to write down the Bethe equations it is useful to define spectral parameters 
x^,u^,u which are interrelated by 

x{u) = ^ + , = xiu^) , M± = M ± i , (2.5) 

and the coupling constant g is proportional to the square root of the 't Hooft coupling 
constant 

These parameters are assigned double indices, i.e. the Bethe equations are written in terms 
of xff^, u^f^ where u = 1, . . . , 7 denotes the nesting level (or node- number of the Dynkin 
diagram), and k = 1, . . . , Ky labels the Bethe roots of level v. 

The Bethe equations for the specific excitation pattern (12.21) are given by (c/ [8]) 

2L-2 + 



.J 



2L-2 

A U \ T — r ^ A I, 



~ I K- 1 _.2/^- ^+ ^ (»4,fc,»4,,) 



1 5'^/^4^A;^4,j 



•^4,fe / ^4,fc •^4J y / •^4,fc-^4,j 

2L-3 - L-1 - L-3 



fj[ a^Jfc - X3,j y x+fc - f}^ 1 - ^V<fca;7,. 



_L-i - _ L-3 1 _ 2 / - 

^3,j T-r I, Xsj -p-j- i g /X^ uXtj 



L-2 i 2L-2 _ + 

^ TT '^5,k ~ UQ j + 2 -rr X^^k X^j 

jj[ «5,fc - - I X^^k - X^j 

L—2 L~l ' L—3 

^ _ TT UQ^k — Uqj — i T-r Ue^k — + t TT '"6,fc — + I 

W6,fc - Ue,j + i JJj U6,k - U5J - \ JJ[ U6,fc - M7,i " f 
L-2 , j 2L-2 1 2 / + 

where the dressing factor a'^{uk,Uj) = e^'>-^'<^k,uj) jg gjygj^ p xp _ 

The (asymptotic) anomalous dimensions of an operator can then be deduced from the 
energy E{g) of the long-range spin chain by 

A = Ao + 2/E((7), (2.8) 

where Aq denotes the operator's classical dimension. In the case of the field strength 
operators (11.61) it is given by Aq = 2L. In turn, this energy is found from the "momentum- 
carrying" Bethe roots living on the central (i.e. 4-th) node as 

2L— 2 

E(9)-Y.\ir--^\- (2.9) 



k=l 



^4,fe 



Let us note that the state JF^, which is a tensor product of 

1^) = bIblcHcHiO) , (2.10) 

has mirror- inverted excitation numbers and hence the same Bethe equations fl2.7p . 

We will firstly solve the one loop-problem at finite L, and secondly construct the ther- 
modynamic limit L oo. Subsequently we will focus on the all- loop equations and derive 
an integral equation for the corresponding root density. 



3 The One-Loop Problem 

The one-loop Bethe equations for the field strength operators Tr JF^ are obtained from the 
asymptotic all-loop equations (12. 7p by taking the hmit (? ^ 0. This yields 

2L-2 i 

7=1 ^3,fc - U^^j + f 
I i \ L 2,L — 2 . 2L — 3 j L — 1 i 



7,, , _ . _ 11 _L i 11 



M4,fc - 2 / JLi ^4,fc - M4,j - ^ ^J^ ti4,fe - M3,j + 2 ~=i "4,fc " ^iSJ + 2 
^ 2 ■ 2Z/ 2 



j^]^ ^5,fc — ^6,j — 2 j=i ""S.fc — M4J -I- 2 
Z/— 2 L— 1 ■ Z/— 3 

^ _ TT Me.fe ~ ""ej ~ TT ""^.fc — ^sj + t TT ""e.fc — mtj + f 

^^TT ^7..-n6. + t _ (3.1) 

7=1 "7,fc - M6j - t 

These can be further simplified with a suitable dualization, see e.g. |32] and references 
therein. We start by introducing the polynomial P[u) for the Uj, roots as 

2L-2 2L-2 2L-3 

p(«) = n - ^4,. + f ) - n - "4,. - f ) = m - 2) n (« - «3,,) . 

It is straightforward to find 

'2^j^—'2^ 2Z/ 3 

P{uA,k + |) _ TT '^■i,k - U4J + ^ _ TT '^i,k " U^j + | 

Piu4,k - I) M4,fc - M4J - i 11 M4_fc - - | ' 
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Applying the same procedure to the roots Uj and Uq we reach an effective system of 
equations, namely 

M4,fe - f / M4,fc - M5,i + I 

L-1 , . 2L-2 i 

By introducing the polynomial Qa{u) for the momentum carrying roots U4 as 

L-l L-l 

Q,{u) ={u+ - "5,. +!)-(«- i)^ - "5,. - 1) 

2L-2 

= ^(2L-1) J](M-n4,,) (3.4) 

i=i 

and noting that 

Qi{u^,k + f) _ / ""S.A: + A TT M5,fc " ^^5^ + ^ _ TT ""S.fc " M4,i + f 

we find that fl3.3p simplifies enormously to the free equation 

1. (3.5) 



U5,k - i 

This is solved by 

M5,fc = ~^ ) ' A; = 1, . . . , L - 1 . 
Plugging this solution back into equation (13.21) we infer that ^4^^ are roots of the equation 

u + ^^M-COt^^' 



n ^ . (3.6) 



"=>-cot(^) + t 
On the other hand, the roots are zeros of the polynomial Q^{u) 

Q,{u) = {u + i)'^-{u-if, (3.7) 
such that we can rewrite (13. 4p as 

Q^{u) = {u+ '-f Q,{u '-f Q,{u - . (3.8) 

We have therewith derived a polynomial Q4{u) of degree 2L — 2 

Q^iu) = {u + i)^ {{u + |«)^ -{u- l^)^) + (m - f )^ {{u - - (m + 1^)^) , (3.9) 

whose zeros correspond to the exact one-loop roots U4 for arbitrary, finite L. All roots of 
this polynomial are real (for a proof, see Appendix 1X1). 
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3.1 Energy and Higher Conserved Charges 

Since the u^^k are the roots of (13.91) we can write 

2L-2 

Q^iu) = -2L(2L - 1) JJ (m - U4,j) . 

i=i 

This allows us to express the energy E-fL{g) = Ejtl + 0{g'^) (I2.9p in the one-loop approxi- 
mation EjrL in terms of Qi{u) 

On the other hand, using the explicit form (13. 9p we derive 

Ql(t) 3.^ QU-t) 3.^ 

-4- = — iL, — ^ = -zL, 

Q4(|) 2 Q4(-t) 2 

and thus verify that the anomalous dimensions (12.81) of the operators TrjF''" are indeed, to 
one loop [291 E], 

A^L = 2L + 6L/ + 0(c/^). 

Incidentally, it is straightforward to generalize the above method to all higher conserved 
spin chain charges, by noting that 

(r- l)(r-2)! (^l^g^^^L^^^ - ^ logQ4HL^-.) 

= ^7^(l + ^)(^^^^-Hr^)- (3-11) 

It is clear that all odd charges vanish. Note that the in (3.20) of [32] for the "Beast" 
diagram are given by ^. 

3.2 Generating Polynomials 

Starting from the set of all-loop Bethe ansatz equations (12.71) and dualizing the and 
roots [8] we obtain 

2 

' + \ -'^ 2L-2 - +1 ^ 2L-4 + 

^ = 11 1 K.,«4,,) 11 _ (3.12) 

^ 2 _ ^ 2Z; 2 ~ -|- 

l = TJ "^'^ ~ ^ TT ^^'^ ~ ^^'^ (3 13) 

j}^ U5,k - Uqj - i Al ^^^^ _ 3™-^, 
L-2 . 2L-4 ~ j 
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Here, Ms are the 2L — 4 roots dual to u^. Note that they are not given by fl3.7p anymore. 
It is however an easy task to derive the polynomials which generate the one-loop roots for 
this set of equations. One of them, generating the M4 roots, is already known, see (13.91) . 
The polynomials generating the and the uq roots are given, respectively, by 

Q5iv) = 3v^^ + {-i + v)^{-2i + v)^ + {2i + v)^ {{t + v)^ + {-2t + v)^) 

-v^ {{-i + v)^ + + v)^ + 2{-2i + v)^ + 2{2i + v)^) (3.15) 

and 

Qeiw) = {w + |z)^ + iw- |i)^ -iw + \%f - (w - \%f . (3.16) 
Some properties of these polynomials are studied in Appendix |X1 



4 Thermodynamic Limit 

In this section we will demonstrate how to construct the thermodynamic limit and how to 
find the root densities using the generating polynomials defined in (13. 9p . (13. 151) . (13. 161) . We 
will present the method, taking Q%{w) in (I3.16P as an example. 
Let us define the density of roots as 

1 

/'i,6(^) = 7$^'5(^-w.)- (4.1) 



L 

i=l 



Using this definition we can write (I3.16P as 



Qeiw) = ]^(w - t^i) = exp J pLfi{z) \og{w - z)dzj (4.2) 



i=l 

and easily obtain an integral equation for the density 



dz = - = Vl{w) . 4.3 



Taking the limit L —>■ 00 and using the identity 

1 1 



X + ie X — ie 

we derive the formula 



-2niSix) , (4.4) 



where 



P6(^) = ^ (Viw + te) - V{w - te)) , (4.5) 



(m) = lim pLfi{u) and V{w) = lim Vl{w) . (4.6) 

L— >oo ' L— >oo 
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It follows from the above formula that V{w) is necessarily discontinuous across the real 
axis (i.e. it has an infinite cut there). To find V{w) we compute 



_ {w + + {w- 10^-^ -{w + 1^)^-^ -{w- 1^)^-^ 

{w + + {w- -{w + -{w- 

Taking the limit is straightforward 

, lm(w)>0 



(4.7) 



V{w) = lim Vl{w) = <{ ^ ' ^ ] ; ^ . {U 
Using (14. 5 p we find for the density 

1 3 

Peiw) = lim PLfiiw) = . (4.9) 

L^oo ZTT + I 

The same procedure can be used to find the density for the momentum-carrying roots 

P4(w) = 7^ (-j^ + , (4.10) 



which is a continuum limit of 

2L-2 



PL,4(^) = 7 5^ 5(^ - «.) . (4.11) 



L 

i=l 



As a validity check we can show that 



oo 



+ (^'-^^ - nr^) > (4-12) 



r 



which is in agreement with the earlier result (13. lip . 

Let us now direct our attention to the density of roots of flavor u^. From the analysis 
of the generating polynomial we infer that these roots form two strings along ±| (see 
Appendix [B] for details) on the complex plane and we have to integrate over two contours 



dv+ + y v~ = A{u), (4.13) 



where A(m) is given by 



"T^ ~\~ — To^ 1 Im(u)> ^ 



A{u) = lim { dulogQ^^L^u) = < ^ + ^ , < Mu)< i . (4.14) 



11 



The solution of (14.131) is 

p,{v±^)=p,{v). (4.15) 

Hence, two roots form a stack with one root of type uq in the center such that 
= (m6 ± |). The concept of stacks was introduced and studied in the context of the 
nested Bethe ansatz of "ferromagnetic" root distributions [32]. Here we find the same 
stack picture in the context of a long operator akin to an "antiferromagnetic" state. The 
emergence of the stack picture for our case is rigorously established in Appendix [Bl 



4.1 Asymptotic All-Loop Effective Bethe Equations 

In the thermodynamic limit only the centers of the stacks receive quantum corrections, 
i-e. u^,k{g) ~ Ufs,k{g) + f and u^,k^-2{g) ~ u^,k{g) - V2 for /c = 1, L - 2. The effective 
set of Bethe equations then readcl 



■^4,fc / j=i ■^A,k -^AJ ^ g I ■^A,k-^A,j 



X 



L-2 + _ - _ 2 1- 



TT ■^A,k -^GJ y I •^A,k-^6,j (4 16) 



where we have used the identity 



L-2 . 2L-2 - _+ i_2/ - - 

■J-|- U6,k — Uqj + I ^6,A: ^4,j 5' /^6,k^A,j ^^s^ 



y^k ~ ) - 1 „2 /™±^± 



1 - g^/x,^ X 



Note that (13.131) splits up into two equations with k running from k = 1, . . . , L — 2 and 
k = L — 2 + l,...,2L — 2 respectively, which can by multiplied with each other to get 
( I4.17p . while equation ( I3.14p is identically satisfied. These effective equations describe the 
scattering of the roots M4 with the centers of the stacks. Rewriting (14. 161) and (I4.17P as 

^tk\ '^TT UA,k - UiJ -i A ^ g"^ / ^A,k^A,j\ 2/ A 

,^4,fc / fJl «4,fc - U4,j +1 \l- gVx^f^xlj J 

UA,k - uej -il- gyx+j^x^j 1 - gyx+^x+j 



"'The use of instead of "=" indicates that these equations are strictly speaking not exactly valid 
for finite L due to differences of 0{j-) for the generating polynomials. See Appendix [Bl for more details. 
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L-2 . 

^\ ue,k - Uqj - i 



2L-2 

}J[ W6,fc - M4,i + i ' 1 - ^7 a;6,fcX+^. 1 - gy x+,^xtj ' 
and taking the logarithm of these equations leads to 

2L-2 



- (logx+fc - logX4 J ^ 27r— + If^*^ ^ 




+ T7 > log ^ ^ : + log ; 1' + log ; , (4.20) 

where the dressing phase shift 2 z ^ is given in (11.11) , and to 

M6,fc - Me.i + i 




The ambiguity in the choice of branch of the logarithm is encoded in the mode numbers 
rik and m^, for which we introduce mode functions = n{uii^k) and = m{uQ^k)- 
As usual we will now replace the sums by integrals as L ^ oo and introduce new variables 
^ _ n(u4) ^ _ m(ue) ^ rjj_^gy g^YYow US to define the all- loop excitation densities through 

p{u) = —^^^^7^ and ri(w) = — One can show that this is consistent as the correct one- 
loop densities can be derived from the effective equations. Finally we take the derivatives 
w.r.t. u and w, respectively. Thus, (l4.2Up and (14.211) become: 

' — 2 TT p[U) 




Ag^ ^ul - Ag^ J J_^iu-u'y + l 

+2. du' Pin') - (lo, +^oiu, 



+2 



r]{w)dw 
{u — wY + 1 



00 

00 ^ / 1 .2 



/ , , 1 — q/x iu^x^iw) , 1 — q/x (u)x (w)\ 

/ dwr]{w)— (log V ^ _ +iog Y/ + + 4.22 
du \ l — g^ix+[u)x [w) 1 — g''/x+[u)x+[w) / 
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and 

rj{w')dw' p{u)du 



= 2 7171(10) + 2 / , \, 2 / 



[w — uV + 1 



oo 



I , , ^d \ — q^lx^{vS)x (u) , l — q^/x {w)x {u)\ 

4.2 All-Loop Energy 

Our strategy is now to solve fl4.23p for the auxiliary root density rj as a functional of 
the energy-momentum-carrying main root density p. Substitution of this solution into 
f l4.23p then yields a single, closed equation for the latter. This is easily done by Fourier 
transformation techniques, see in particular [T5] . 

The Fourier transforms of the densities p{u) and ri{w) are defined as 

/oo roo 
e''''p{u)du, f]{t) = e-^ I e''"" r]{w)dw , (4.24) 
-oo J — -oo 

where we have included a factor e"'*'/^ for convenience. Including the same prefactors, the 
double-Fourier transforms of the kernels 

K{w, u) = — d^ log _ , , 4.25 

2711 1 — g^/x [w)x+{u) 

1 , 1 — / x~ (w)x~ (u) 
i/(u^„) = -a„log^-^5y-i-Ai. (4.26) 

may be explicitly computed as, respectively, 

KM = 2.9^(1 -sign ,,)|,|,H.IHn^°(2gW)^.(2^J^'1)-_^o(^^ 

Hit, f) = - 2./(l + sign ,,)|,|,H.|Hn^o(2gl'l)^.(2gl';i) + ^»(2gl«1)^.(2gM) . (4.28) 

^ ' ' ^ ^ ^ ^' ' 2g{\t\ + \t'\) ^ ' 

Fourier transforming (14.231) diagonalizes all terms containing r], and one obtains 
= e\'\ flit) + flit)- pit) 

/ dt' p{-t')e\'\k{t,t')e\''\ / dt' p{-t')e\'\H{t,t')e\''\ . (4.29) 

As announced above, we may now solve for fiit) in terms of an integral transform of p{t) 

^(^) = ;r^('^W-4^'^/ dt' k,i2gt,2gt') pit')^ , (4.30) 

whose kernel KQit,t') is defined as in (11.31) . 
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Now we treat fl4.22p in a likewise fashion. The Fourier transform of the Lh.s. is given 

by 

I r°° , f 1 1 



, du , , e**" = e-l*l/Vo(2^t). (4.31) 

Thus, equation (14.221) reads in Fourier space 
Jo{2gt)=e'p{t)-p{t)+fi{t) 

+^9^t dt' {k^{2gt, 2gt') + ka{2gt, 2gt')^ p{t') 

POO 

~AgH / dt'k^i2gt,2gt')fiit'), (4.32) 
where we have used (11.11) , and Ki is given in (ll.4p while Km = Kq + Ki is 

A-„(M-)^''''""°"'>:;f°'" (4.33) 
Now we may substitute (I4.30p into (I4.32p . This leads to a single equation for p{t) 

p{t) = e-^\l + e')(^Jo{2gt) 

+AgH / dt'(^—Koi2gt,2gt') + ^—TKii2gt,2gt'))p{t') 
Jo ve+i e+i / 

-4gH ^ dt' {Km{2gt, 2gf) + k^{2gt, 2gf) + kd{2gt, 2gt')^ p{t')^ , 

(4.34) 

as announced in (II. 8p at the beginning of this paper. Since we have two degrees of freedom 
per unit length, i.e. the state is doubly-filled, the density (I4.34p is correspondingly normal- 
ized to 2, i.e. /3(0) = 2, see (I4.24p . The dressing kernel Kd{t,t') is written in (II. 2p . and we 
just proved that the nesting kernel Knit, t') is indeed given by the structurally very similar 
expression (11.90 . In fact, 2Kn and K^ differ only by a single sign in the denominator of 
the diagonal kernel t/(e* ± 1). 

The energy of the state can now be computed from the density by Fourier transforming 
( IXQj) . It reads 

E{g) = AL dt pit) . (4.35) 

This leads via (12. 8p to the anomalous dimension (ll.7p announced at the beginning. 

We expect our solution to be valid at arbitrary values of the coupling constant. The 
integral equation (I4.34p is however (presumably) too complicated to be explicitly solvable. 
Being of Fredholm-type, it is however very easy to expand in small g to high orders. We 
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thus find for tlie perturbative energy to e.g. seven-loop order 
L 4^ 4^ 164 '-vyiy 



/632661 1665 ^,^, 8 



/29056443 87525 , , , , s i m 

+ — — C(3) + 8505 C(5) + 5670 C(7) g'' + 



256 



/351914817 2244573 , ,2 



+114723 C(5) + 90909 C(7) + 63504 C(9) J g'^ + . . . . (4.36) 

We notice tfiat zeta functions of odd, but not even, argument enter the energy. This is not 
surprising, as these are directly generated by the dressing kernel fll.ll) . c/jH]. If we assign 
a "degree of transcendentality" k to C{k) we see that the contributions at a given loop 
order are, in contradistinction to the case of large-spin twist operators, see [IHl [lH [I9] , 
not of constant degree. One may nevertheless observe that at a given order / the degree is 
bounded, and always saturated, by 2/ — 5. Note also that all zeta- function coefficients, as 
well as all rational numbers, turn out to be integers after factoring out inverse powers of 2. 

Conversely, the dressing factor is the only source of (^-function terms in the expansion 
f l4.36p . By this we mean that dropping the dressing factor from the asymptotic Bethe 
ansatz equations (12. 7p would eradicate all terms containing (^-functions in (I4.36p . and would 
generate only rational loop contributions. The latter would precisely agree with all terms 
of "transcendentality degree zero" in (I4.36p . 

It is also interesting to note that the thermodynamic expansion of (I4.36P does not 
coincide with the energies of finite length operators (see Appendix [C]) , as opposed to the 
case of the pseudo- vacuum state TrW^ [211 [33]. By this we mean that the exact anomalous 
dimension of the operator Tr JF^, even below wrapping order, is not exactly proportional to 
L, while it is for TrU^. This is very likely due to the length changing processes starting for 
field strength operators at two-loop order, i.e. the number of fields in a local operator is not 
a conserved quantity at higher loops. Correspondingly, and in contradistinction to a BPS- 
state Tr Z or a fermionic pseudo- vacuum TrW^, the one-loop field strength operator Tr J-" 
is not an exact eigenstate and will pick up higher-order quantum corrections: Tr J-'^+0{g^). 



5 Maximal Filling and Nesting 

We would like to demonstrate that the emergence of a nesting kernel (II. 9p . which we 
argue to be closely analogous to a dressing kernel (ll.2p . is a rather generic mechanism 
if two prerequisites are met. The first is that the states satisfy a special maximal filling 
condition, and the second is that the state is irreducibljl^ nested. 

^ By this we mean that the nesting cannot be removed by an exact duahzation of the Bcthe roots. 
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Let us first study two simple cases where only the first, but not the second prerequisite 
is fulfilled. These are the "highest energy" antiferromagnetic state of the su(2) sector 
[26l and the fermionic pseudo- vacuum state TilA^ [HI], [33]. As a by-product we will 
find the correct perturbative expansion of these states, as the expressions in [26l [311 [31] 
were obtained with a trivial dressing factor cr^ = 1, and need to be revised starting at four 
loops. 

For the su(2) antiferromagnet the occupation numbers for an even length L operator 

are 

(i^i, K2, Ks, K,, K,, K,, Kj) = (0, 0, 0, |, 0, 0, 0). (5.1) 

It is straightforward to repeat the analysis of [261 presence of the non-trivial 

dressing phase (11.11) . and one finds that the Fourier-transformed density of roots satisfies 
in the thermodynamic limit the linear integral equation 

m = ^^(^M29t)-'^9't dt' ka{2gt,2gt') pit')^ . (5.2) 

This is a half-filled state, which is the maximally possible filling in this sector. Correspond- 
ingly, the density is normalized to 1/2, i.e. /)(0) = |, see fl4.24p . The energy is as always 
given by fl4.35p . Dropping the convolution term on the r.h.s. of (15. 2p yields the density 
for the Lieb-Wu ground-state energy of the fermionic Hubbard model in a closed form, 
see [26]. The 5u(2) sector of A/" = 4 gauge theory is a supersymmetric deformation of the 
latter. It is described by adding the backreacting convolution. Apparently, the density 
pit) can no longer be found in closed form. However, it is interesting to work out the weak 
coupling expansion of our model from (I4.35p .( !5^ . One finds to the first few orders 

^ = 2 Ca(l) - 6 + 40 Ca(5) - (^350 Ca(7) + 32 C(l) Ca(3) C(3)) / 

+ (^3528 Ca(9) + 96 W^f C(3) + 320 Ca(l) Ca(5) C(3) + 320 Ca(l) Ca(3) C(5)) g"" 

38808 Ca(ll) + 1600 Ca(3) Ca(5) C(3) + 3360 Ca(l) Ca(7) C(3) 

+ 1024 Ca(3)' C(5) + 3264 Ca(l) Ca(5) C(5) + 3360 Ca(l) Ca(3) C(7))^ 

+ (^453024 Ca(13) + 6400 U^f C(3) + 15680 Ca(3) Ca(7) C(3) 

+37632 Ca(l) Ca(9) C(3) + 512 Ca(l) Ca(3)' C(3)' + 17792 Ca(3) Ca(5) C(5) 
+34944 Ca(l) Ca(7) C(5) + 11200 W^f C(7) + 34944 Ca(l) Ca(5) C(7) 

+37632Ca(l)Ca(3)C(9)V'' + -- - • (5-3) 
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Here, we have decided to distinguish the alternating "fermionic" Ca-function 



n=l 

from the ordinary "bosonic" ^-function 

r(A;) Jq t e* - 1 f^^n'' 

even though they are related by the formula Ca{k) = (1 — 2^^'') ({k). Note also that 
^a(l) = log(2). Using these relations we can "simplify" the expression (15. 3p to 

^ = 2 log(2) - I C(3)/ + y C(5)^^ - (24 log(2) C(3)^ + C(7)) / 

+ [54 C(3)3 + 540 log(2) C(3) C(5) + C(9)) / 



/ 11655 
( 1701 C(3)^ C(5) + 3060 log(2) C(5)'^ + ^- log(2) C(3) C(7) 

128 V 



(71 505 
288 log(2) C(3)^ + 18135 C(3) C(5)^ + ^ C(3)^ C(r) + 67158 log(2) C(5) C(7) 

57Q79Q1 5 \ 

+65709 log(2) C(3) C(9) + C^jg'' + ■■■■ (5-6) 

This however obscures the distinction between the contributions stemming from the dress- 
ing factor (the ^-terms) and from the fermionic Hubbard model (the Ca-terms). In fact, 
omitting all terms containing ( (in (I5.3p . but not in (I5.6p !) leads back to the Hubbard 
ground-state energy. This admixture of bosonic ( and fermionic (a is further evidence 
that planar AdS/CFT is a supersymmetric generalization of the purely fermionic Hub- 
bard model as employed in [20] ■ Note also that the coefficients multiplying the (- and 
Ca-functions in (15.31) are all integers. Interestingly, we see that the terms in (15. 3p are still of 
constant degree of transcendentality (2/ — 1) at a given loop order / if we assign a "degree 
of transcendentality" k to both ({k) and Ca{k)- 

We now turn to the su(l|l) state TtU^. The occupation numbers are 

{K,, K2, K3, K,, K5, Ke, Kr) = (0, 0, 0, L, L - 1, 0, 0). (5.7) 

In the picture of Fig. 1, we first replace all L fields Z in the BPS vacuum Tr by L bosons 
X, and then turn the bosons into fermions U. This state may be dualized, and the two-level 
nested Bethe equations may be converted to a single level [H]. Extending the analysis of 
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to the case of a non-trivial dressing phase (II .ip . we find that the Fourier-transformed 
density of roots satisfies in the hmit L — > oo the equation 



m = e-* (^Ulgt) -'^9^t dt' {k^{2gt, 2gt') + 2Ka{2gt, 2gt'))p{t'\ 



{5.1 



This is a filled state, with one excitation per lattice site, which is the maximally possible 
filling in this sector. Correspondingly, the density is normalized to 1, i.e. p{0) = 1. Working 
out the weak coupling expansion of this state from fl4.35l) . fl5.8l) . one finds to e.g. eight-loop 
order 

^M. = 2 - + 58^^ - (^518 + 32 C(3)^ / + (^5228 + 480 ((3) + 320 C(5)^ / 
57280 + 6144 C(3) + 4928 C(5) + 3360 C(7)^ c/^° 
+ ( 665344 + 76768 C(3) + 512 ((3)^ + 64960 C(5) 
+52864 C(7) + 37632 C(9)^ g^^ 

8070352 + 965856 C(3) + 13312 ((3)^ + 833792 C(5) + 10240 C(3) C(5) 
+713056 C(7) + 602112 C(9) + 443520 C(ll) )g^^ + ... . (5.9) 



The (^-functions are exclusively generated by the dressing factor. We note similar "tran- 
scendentality properties" as in the Tr J^^ case of the last section, namely that at a given 
loop order / combinations of zeta functions with odd arguments occur up to and including 
degree of transcendentality 21 — 5. Again all zeta-function coefficients, as well as all rational 
numbers, turn out to be integers. 

Let us finally present a third example, namely a certain 5o(6) singlet state. At one- loop 
this state is the highest energy state of a 50 (6) magnet [H |35]. At higher loops, however, 
the 50 (6) subsector is not closed anymore and thus one is forced to use the full p5u(2, 2|4) 
Bethe equations. The excitation scheme for this state reads 

{K,, K2, K,, K,, K,, K,, Kj) = (| - 2, ^ - 1, |, L, |, | - 1, - 2). (5.10) 
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After the dualization of the and M3 roots one is left with the following equations 

, ■L-2 ~ i 

f-1 ~ , i L ~ + 

4,fe \ -TT 4,fc ■^4,j ^4,fc^4,i 2/ \ TT 4,fe "^Sj -pr -^4,^ "''Sj 

— =11^ —- {Ui,k,U4,j) \\— ^11^ — 

Z ^4,fc''4,i 

U5,k - U6,j + t -TT ^5,fe - 2:4 



. _i n 



'^5,k — — ^ X5^k — X^^j 

L 
2 



^ 1 , .L-2 



For the highest energy state we may assume M2j = u^j and uj,^k = U5,k for j = 1, ■2 — 1 
and k = 1, L—2, respectively. After this is done one notes a striking structural similarity 
to fl3.12l) - fl3.14l) . Apart from the different overall number of magnons the resulting equations 
differ only by the power of the interaction term in the equation for the M4 roots. It is thus 
plausible to assume that the M5 and uq roots will again form stacks, even though this seems 
to be more difficult to prove in this case. Under this assumption, the derivation of the 
integral equation for the principal density is straightforward and follows the same lines as 
above. In Fourier space one now gets, in great similarity to fll.Sp . 

-AgH dt' {-^K^{2gt,2gt') + ^—^K^{2gt,2gt')\ p{t') 

-^9^i dt' (2kr,{2gt,2gt')+kd{2gt,2gt')^ p{t')^ . (5.12) 

Interestingly, the very same nesting kernel Kn ( \1.9\\ appears, with an overall factor of two, 
as in the case of the field-strength operator, cf (11.81) . This is precisely what one should 
expect. 

It is again rather straightforward to find the weak-coupling expansion of the energy to. 
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say, four loops (we have highlighted the terms generated by the dressing phase) 

^ = 2/3(1) + C„(1) 

4/3(1) m + 6/3(3) + 2/3(2) (»(!) - /?(1) Ca(2) - ^ Ca(l) Ca(2) + ^ Ca(3))y' 

+ (8(3{1) (3{2f + 8/3(2) /3(3) + 24/3(1) /3(4) + 40/3(5) + 4/3(2)^ ^(1) 
+12/3(4) C„(l) - 4/3(1) /3(2) C„(2) - 2/3(3) Ca(2) - 2/3(2) Ca(l) Ca(2) 

+ i /3(1) Ca(2)^ + \ Ca(l) Ca(2)^ + /3(2) C„(3) - \ Ca(2) Ca(3) - ^ /3(l) Ca(4) 
-7Ca(l)Ca(4) + 7Ca(5)V' 



- (^16/3(1) /3(2)3 + 16/3(2)2/3(3) + 16/3(1) /3(3)2 + 96/3(1) /3(2) /3(4) + 44/3(3) /3(4) 

+40/3(2) /3(5) + 200/3(1) /3(6) + 350/3(7) + 32/3(1) /3(3) C(3) + 8/3(2)^ Ca(l) 
+8/3(3)2 ^^(1) ^ 48^(2) /3(4) Ca(l) + 100/3(6) Ca(l) + 16/3(3) Ca(l) C(3) 
-12 /3(1) /3(2)2 C„(2) -8/3(2) /3(3) Ca(2) - 24 /3(1) /3(4) Ca(2) - 10 /3(5) Ca(2) 
-6 /3(2)2 C„(l) Ca(2) -12/3(4) Ca(l) Ca(2) + 3/3(1) /3(2) Ca(2)2 ^ ^(3) ^^(3)^ 

+ 1 m) Ca(l) Ca(2)^ - \ /3(1) Ca(2)^ - ^ Ca(l) W^f + 2 /3(2)2 (,(3) 

+4 /3(1) /3(3) Ca(3) + y /9(4) Ca(3) +4/3(1) Ca(3) C(3) + 2 /3(3) Ca(l) Ca(3) 

+2 Ca(l) Ca(3) C(3) - m Ca(2) Ca(3) + ^ W^f Ca(3) + | /3(1) Ca(3)' 

+\ Ca(l) Ca(3)' -6/3(1) /3(2) Ca(4) - ^ /3(3) Ca(4) - 3 /3(2) Ca(l) Ca(4) 

+ ^ /5(1) Ca(2) Ca(4) + \ Ca(l) Ca(2) Ca(4) - ^ Ca(3) Ca(4) + \ /3(2) Ca(5) 

Ca(2) Ca(5) - y /3(1) Ca(6) - ^ Ca(l) Ca(6) + Ca(7) + ■ ■ ■ ■ (5.13) 

For this state, an interesting new set of numbers appears, namely the Dirichlet /3-function 
evaluated at positive integers: 

If k is odd this leads to vr*"' times rational numbers (related to Euler numbers). If k is 
even the numbers /3(/c) cannot be expressed through tt's. /3(2) is Catalan's constant. All 
coefficients multiplying the products of C,-, Co, and /3-functions are integers after factoring 
out inverse powers of 2. 
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Note that fl5.12p is the higher loop generahzation of a one-loop result worked out in [T] . 
The latter is of course reproduced from (15.131) since (3{1) = 7r/4 and Ca(l) = log(2). Note 
that this state's energy satisfies, just as the half- filled su(2) state's energy (15.31) . a constant 
transcendentality principle. 

6 Outlook 

We have shown in detail how to compute the anomalous dimension of the field strength 
pseudo-vacuum (II. 6p in the thermodynamic limit from the nested asymptotic Bethe equa- 
tions of [8]. Several techniques to reduce the number of equations were introduced and 
a single effective integral equation (II. 8p for the distribution density of Bethe roots was 
derived from a starting set of five. Combining this equation with the expression (11.70 re- 
lating the density to the operator dimension, it is straightforward to find the weak coupling 
expansion of the latter to any desired order, cf (14.360 . Incidentally, as our equation is an- 
al5d;ic in the vicinity of g = 0, it should also be, by analytic continuation, just as valid at 
any value of the coupling. It would be very interesting to analyze it in the strong coupling 
limit g ^ oo, and to interpret the corresponding state in string theory. The techniques 
developed in [36] might be useful here. 

Interestingly, the influence of the nesting on the effective Bethe equation results in a 
kernel (11.90 which strongly resembles the dressing kernel (11.20 recently proposed in [Ti] . 
This leads us to suggest that the dressing phase should originate from the elimination of 
further, yet to be found auxiliary Bethe roots. A formal procedure, unfortunately plagued 
with difficulties, is briefly discussed in Appendix |D] (see also The detailed mechanism 

for how this happens therefore remains to be worked out. The techniques presented in the 
present paper might prove helpful in this respect. 

We would also like to stress that "extra Bethe ansatz levels" in order to improve the 
asymptotic equations of [S] have been proposed previously. These are natural from general 
arguments concerning finite size effects [3H], from concrete indications that the asymp- 
totic Bethe ansatz needs to be corrected [39], and finally from studies indicating that the 
BDS/Hubbard magnon dispersion law at strong coupling PU] is to be corrected in a finite 
volume [H] (see also [12] and [13]). In particular, Hubbard-type models are able to create 
long-range integrable systems from short range interactions [2S] (see [H] for a detailed 
study on the relation of the nested to the effective Bethe equations near the antiferromag- 
netic vacuum). Furthermore, similar mechanisms have also been proposed on the level of 
the string sigma model [271 EH] • It is interesting to note that one key feature of [26l [271 EH] 
is that the BPS vacuum corresponds to a non-trivial distribution of top-level Bethe roots. 
This qualitatively agrees with our result, which suggests that the "physical" BPS vacuum 
states (I1.5P are created by filling up a truly empty "unphysical" reference state. See also 
the closely related comments in [l5]. 

The number of hidden Bethe roots creating the dressing factor of [11] is infinite if the 
ansatz is to be asymptotically exact for short operators. This should be related to the 
non-compact nature of the AdS/CFT system. In such a situation Bethe equations might 
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not necessarily furnish the most effective description, and an approach based on Sklyanin's 
separation-of-variables technique might be more appropriate. This would then replace the 
Bethe equations by functional equations for a nested set of Baxter-Q functions. Some of 
the latter should be non-polynomial in nature. Possibly the techniques used in ^61 HTj 
might be useful in this context. 
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A Properties of the Generating Polynomials 

In this section we will analyze some properties of the generating polynomials Q4 and Qq. 
The mathematical description of the formation of stacks will be given in the next section. 

Let us start with the following observation: if Q{x) is a polynomial with only real roots 
then 

S{x) = Q{x + i s) + a Q{x — is) 
with |a| = 1 also has only real roots. To prove this one observes, that for any real y 

n 

\Q{x + iy)\ = \Qix - iy)\ =Y[\ix - Wj) + iy\ 

i=l 

grows monotonically with increasing \y\. Furthermore, if Xj satisfies: 

S{xj) = Q{xj + i s) + a Q{xj — i s) = 

then \Q{xj + i s)\ = \Q{xj — i s)\. Together with the previous observation this implies 
Xj G M. Applying this theorem twice to one proves immediately that Qq{w) has only 
real roots. Similarly 

Q^{u) = W{u + ^) - W{u - '-), Wiu) = + z)^ -{u- ^)^) 

also has only real roots. If one orders the roots of Q{x) as follows xi < X2 < ... < Xn one 
finds that 

Q'{xi) g'(x,+i) <o. 

Thus, either Xi or Xj+i are zeros of Q'{x) or there exist a 6 with < < 1 such that 

Q'{xi + 6{xi+i - Xi)) = 0, i = 1, n - 1 . 
In case of Qg{w) one finds 

-^Qeiw, L) = L Qeiw, L - 1) . 
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Thus, between any two zeros of Qq{w,L) there hes a zero of Qq{w,L — 1). Furthermore 
it is fairly easy to prove that the resultant (given by the determinant of the corresponding 
Sylvester matrix) of Qq{w,L) and Qq{w,L — 1) is always non-vanishing. Finally, it is 
possible to find an approximate formula for the extreme roots of Qq{w, L). Expecting this 
roots to scale with L we set 



Using 



lim (1 + 4)^ = 



we write 

Q,[±aL, L) = -^{±aLf cos sin + 0{-). 

Setting a = one finds that Q%{^L + e, L) changes sign with e and hence, 

Wmax^±--L, L»l. (A.l) 
TT 

We conclude that for any L > 2 the generating polynomial Qe{w, L) has (L — 2) distinct 
real roots, forming a dense set in R for L ^ oo. Curiously Qq{w,L) is the polynomial 
solution of 

3z 3 2 z z uj 

Qe{w + + Qeiw - y) + Qe{w + -) + Qeiw - -) = 2^ Qe{-) 

This should be interpreted as the Baxter equation for the Uq roots. Similarly one can find 
extreme roots of {u, L) and {v, L). Surprisingly they are also given by ( lA.ip . 



B Formation of Stacks 



Extreme roots of are indeed real (as demonstrated above). In this section we will, 
however, prove that almost all roots of Q5{v, L), for large values of L, occupy two contours 
shifted from the real axis by ±|, respectively. In other words we will show that 

Q,{x, L) = Qe{x+'-,L) Qe{x ~\.L) + 0(i) (B.l) 

where 0{j^) can be neglected in the large L limit. 
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B.l Power Expansion of the Polynomials 

After some computation it is possible to find for the polynomial Q5{v) 

0. (.) = 3.- - 2 ^ (^) co4 - 4 X: 2' (^) cos ^ 

fc=0 ^ ^ k=0 ^ ^ 

+2 (1^''" T 1 + ^ - ^; 2] 

k=0 ^ ^ 

+2 E U^_ y-V^-'cos'^2F^[-L,k-2L■,l + k-L■,2] 

k=L+l ^ ^ 
I.— n V / 



fc=0 

Thus one finds immediately that 

2L-4 



n=0 



where the coefficients are given by 

e„=_2l./ ^eos^P4^-4 2--"l..^ ^ cos " 



2L-n 2 \2L-n 2 



(B.3) 



forn = L, L + 1, 2L - 4 and 



L \_n{2L-n) ^ .^_J L \_n{2L-n) 



c„ = - 2 ( ^ cos - 4 2^^-" I ^ cos 



2L-n 2 \2L-n 2 



.2(^)2-cos^..,-.,-.;l..-n;2].Q^)cos^(^)4-t 

(B.4) 



25 



for n = 0, L — 1. Similarly one finds 

O T 



Qe(w + \)Q,{w - 1) = - 2 f; {^^) w-^-^ cos ^ - 2 X: 2^ (^) cos ^ 

fc=0 ^ ^ fe=0 ^ ' 



+ 2 1; ^^^-^ (^^) cos ^ (^^Fi [-L, 1 + L - ^; 2] 
-2Fi[-L, -A;;l+L-A;;-2]^ 

A:=L+1 ^ ^ ^ 

-2Fi[-L,A;-2L;l + A;-L;-2]^ . (B.5) 
The corresponding coefficients in the power expansion of 

2L-4 
n=0 

can easily be read off 

\ n / 2 \2L — n/ 2 



cos -'IT-^ I 

\n ) 2 V2 

+2 {^^^_ J cos (^2^1 n - 2L; 1 + n - L; 2] 



(B.6) 



-2F1 [-L, n - 2L; 1 + n - L; -2]^ 
for n — L, 2L — 4 and 

\n J 2 \2L-n) 2 

-<./2)-«1¥)-(./2)^^---1¥) 

+2 2^- Q cos {^,F,[-L, -n-l + L-n- 2] 

-(-l)^-^Fi[-L, -n; 1 + L - n; -2]^ (B.7) 
for n = 0, 1, L — 1 . 
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B.2 Emergence of the Stack Picture 

In this subsection we assume n > L {n < L can be analyzed analogously). Furthermore 
without loss of generality we consider L and n to be even. The non-zero coefficients are 
then given by 

Cn{L)^2{-l)^(^^^^_^(^,F,[-L,n-2L;l+n-L;2]^^^ (B.8) 
dn{L) = 2(-l)t (^^_ ] Lfi[-L, n-2L;l + n-L;2] 



-2Fi[-L, n-2L]l + n- L]-2]- 2^^-" 
B.2.1 n = 2L — a case 

The case n — 2L — a with a > 4 and finite a is fairly easy to analyze. Noting that 



(l + L-a)fc ' 'V L 2L2' 

one can immediately shown that for the large L expansion of 2Fi[—L, —a; 1 + L — a; ±2] 

2Fi[-L, -a; 1 + L - a; 2] = 3" (^1 + ^ a{a - 1) + a(a^ + Sa^ - 7a + 3) + ...^ 

2Fi[-L, -a; 1 + L - a; -2] = (-1)" {^-\ «(« - 1) + ;^«(«^ - 7a' + 13a - 7) + ...^ . 
Expanding c„ and (i„ we note that 

for large values of L. 

B.2. 2 Other values of n 

Other values of n are much more difficult to analyze since their generic dependence on L 
is 

n^aL + (5 , l<a <2 

and the summation limit of 2F1 depends now on L. Here we will show that ^ — 1 is at 
least of order 0{j). We start with noting that for this values of n one has 

2Fi[-L, n - 2L; 1 + n - L; 2] > 1 + 2^^-" + ^^'^ ~ '^^] (2 + 2^^-"-^) + ... 

Li — n — i 
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i.e. terms of the form (21-71) " least of all sub-leading. Next we observe that 

L 



2-' > 



j=0 



n f2L 
' =^2ZU 



ah 



n 



This implies that the (^^) term in dBM is at least sub-leading. Let us define the function 

L 



f2L-n,L{^) = 

j=0 



L 

j J \2L -n-j 



= Pi 



(L-ft,-2L-l) 



1 - 2x 



where Pm'^\x) stands for the Jacobi polynomials (see e.g. [IH]) and n = 2L ~ n. It is easy 
to prove the following properties of fn,L{x) 



(B.IO) 



A straightforward application of these is the relation 



which can be used to get an upper bound 

]2F^[~L,n-2L-l + n-L;2] < 2^ 
2L — nJ 



2L 

n 



f2L-n,L{x) is a polynomial of degree 2L — n and thus has 2L — n roots. Because of Descartes 
rule and the properties of Jacobi polynomials all the zeros are real, distinct and negative. 
Furthermore because of (IB.lOp 



L— — 



where all Xj E (—1, 0). We thus can writ^ 



/2L-n,L(2) 



f'. 



2L-ri,L\ 



2Fi[-L,n- 2L;l + n- L;2] 



2F1 [-L, n - 2L, 1 + n - L, -2] 



n 



\x,\+2 



\Xj\ 



\Xj\ 



^For simplicity assume Xj 7^ — i. 
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If the roots are densely distributed the right hand side grows exponentially with growing 
L. To say something more about the root distribution we observe that 



Since Pm'"\y) with a > — 1 has roots densely distributed on (—1,1) we conclude that 
the roots of f2L-n,L{x) are densely distributed everywhere on (— oo, 0) which completes the 
proof. 



C Exact Results for Finite L 

Here we present the energy for the field strength operators fll.6p of finite length L in powers 
of the coupling constant. Their anomalous dimension is fl2.8p . The results were obtained 
using Mathematica. The ^-functions are exclusively generated by the dressing factor. 



L=4: 



L=5: 



L=6: 



L=7: 



^3 = 9 - y + 315g' + ... 



E, = 12- 50/ + 1^,^ - + 216 C(3) 

+ (^^^^ + 3302 C(3) + 2160 C(5) ) / + 



^ 2 ^ 80 ^ V 1600 7^ 

/ 1539949881 8307 C(3) „ x\ « 

+ + ^ + 2700 C 5 + • ■ ■ 

V 32000 2 sv ^ 

837 2 6278355 . /29266837713 . r 

E-g = 18 + - + 324C 3 ]g^ 

^ 11^ 10648 ^ V 5153632 7^ 

/76857234976107 605205 C(3) , x\ s 

+ + ^ + 3240C 5 + ... 

V 1247178944 121 sv 7^ ^ 



179 2 76603 4 /181131695 ^ \ « 

/8959397257 11641 C(3) , 
V 131712 2 sv ; |y 



^Please note that for this particular "short" state wrapping effects compete with the deformation caused 
by the dressing factor. The Bethe equations are not rehable anymore, and we thus omit the 0{g^) term. 
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L=8: 

4380 2 125533809 . /8840715968859 ^, \ ^ 

Eg = 24 / + o^- + 432C 3 0^ 

^ 43 ^ 159014 ^ V 1176067544 7^ 

/346753221469919673 12323484C(3) « 
+ ( 4349097777712 + 1849 +^™<<^T 

There are several conclusions one can draw from these exact results as compared to 
fl4.36p . It is clear that the higher- loop corrections are not exactly proportional to L for 
finite length operators. This is presumably caused by length changing processes, see [8]. 
Also, the transcendentality principle is violated. The terms of highest transcendentality, 
however, are strictly proportional to L and precisely agree with the ones derived in (14.361) . 



D Emulation of the Dressing Phase? 

In this paper we have argued that a nested Bethe ansatz naturally gives rise to a convolution 
structure very similar to the one in the proposed all-loop dressing phase. In this appendixl^ 
we will ask whether we can modify or supplement the asymptotic Bethe equations of ^ 
in order to construct a well-defined set of equations leading to the correct dressing phase. 

We notice that the difference between the nesting phase (11. 9p and the dressing phase 
(II. 2p is just, apart from a factor of two, a minus sign in the diagonal kernel l/(e* ± 1) 
inside the convolution. The factor of two is easily fixed by using two instead of one set of 
auxiliary roots, just as in the case of the so (6) nesting, see (I5.12p . We claim that the sign 
may also be formally fixed by emulating the full set of asymptotic Bethe equations with 
dressing factor by the following set of equations: 

<-yj I \ / T* T' I /'I I IT '"f 



n 



4,fc -^4 J y / "''4,A,-^4,j 



X 



^4,fc / j^i \-^. 

Ki ^ 2 1- K-i - K., - Kr 2 1- 

n-L y /■^4,fe-^l,i TT •^4,k -^Sj TT •^4.,k "''Sj TT ^ il / •^4,fc-^7j 

1 ri'2 I n^i^ ry A A ,7^~l~ ,y J. A ,7^^~ ry A A T 2 / /y»~l~ /y 

^ y /■^4,k-^l,j ■^4,k -^Sj j^i •^4,k -^bj j^i ^ i/ /•^4,fc-^7j 

1 _ n2 /rr~ rp_ Ma - _ Ma - _ Ma 1 _ „2 /„- „_ 

n-L y /•^4,fc-^l,j TT -^4,*; -^SJ TT 4,fc "^SJ TT ^ /•^4,fe-^7j 

1 ri'2 j /yj~l~ ly A A /yj~l~' ry A A /y«~l~ ly A A T 



/"( 2 / /yj"!" rtf .A. .A. /yj /y .X. .X. /y~l~ y .X. .X. T /l2 / /y 



(D.i; 



^ The results of Appendix |D] were known to us for some time, but we did not include them in an 
earlier version of this paper due to the various difficulties discussed below. While preparing an improved 
version of this article similar ideas appeared in [37, . We therefore decided to include our notes in order 
to discuss the problems which would have to be resolved in order to make this and similar mechanisms 
viable. Furthermore, our procedure differs in several ways from |37| . In particular, we distinguish nesting 
from dressing roots (there is no direct interaction between them), and we do not assume a macroscopic 
number of momentum carrying main roots. 



30 




31 



These equations are very similar to the ones in Table 5 of [8j. However, the dressing factor 
Yl cr^ of Table 5 is now missing and is supposed to be emulated by 



Ma 



4,fc y ' l,i 4,A; -^aj **'4,A: -^Sj "^4,*: 9 I -^7. 



7=1 4,fc y /"''Ij -^4,*: -^SJ 4,fc "''Sj 4,A: ^ ' 7,j 



Note that we have introduced an additional set of 2 x 3 x Md auxiliary "dressing" roots 
denoted by a bar on the flavor indices {1, 2, 3}, {5, 6, ?}. They should not be confused with 
the set of 3 + 1 + 3 "nesting" Bethe roots, carrying unbarred indices, with multiplicities 
(i^i, 7^2, -^3, -^4, -^5, -^6; -^t)- Thcsc are identical to the ones in [8j. We will now prove 
our claim, which however requires the Md — > oo limit. We will also show, as required by 
consistency, that the emulated dressing phase indeed factorizes 



as required by the asymptotic Bethe equations of Table 5 in [8]. Notice that there are 
various ways to rewrite the equations fID.ip - fID.Sp since we may always perform dynamic 
transformations [8j and therewith trade back and forth, respectively, roots of type 1, 7 and 
type 3, 5. 

The derivation of the dressing factor formally proceeds as in the main body of this 
article. We will assume the same mechanism of stack formation. We may then eliminate 
the dressing roots of type {1,3} and {5,7} and write "effective" equations coupling the 
dressing roots of type 2 and 6 to the momentum carrying roots of type 4. This yields 



'^fc(^4,l; • • • ; ^i,Kd — W ^ 



-J 1 ~ y^/^4,fc^5^j 1 ~ 9'^ /^i^kX^- 



7=1 ^2,fe - U2J + I J-jL 1 - 9'^|xl^^xX^^ 1 - yV^5;A:<i ' 



-^-^ 1 — xZi,X^ , 1 — Q-^ XjuXj , 



9'^/xlkX-Qj ^-9'^/xXkH,j 



U6,k - UBj + I JJ[ 1 - g'^/x^^xX^ 1 - g'^/xXj^xlj ' 

3^k 



We should however mention that these equations do not have any proper solutions for finite 
values of Md- Nevertheless, one may formally proceed and take Md to infinity. Note that 
we do not assume any of the original quantum numbers (L; Ki, K2, K3, K4, K^, Kq, K-j) of 
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the asymptotic ansatz to be thermodynamically large. The remainder of the derivation 
proceeds again by transforming (ID.lOp into Fourier space 



= e\'\m - i{t) - ^-^ / dt' J2 e-^""^el*l [k{2gt, 2gt') + H{2gt, 2gt')) e'*''/^ (D.ll) 

where k{t,t') and H{t,t') are given by f l4.27p and f l4.28p and i{t) is the density of the 
dressing roots. If we do not assume the distribution of roots to be an even function, we 
end up with two dressing densities, for t > and t < respectively 

a±t) = -^-TTY I ( E e^''''^ kmi^gt, 2gt') + e^^*'"^^^(2^t, 2gt'] 



Y,e^''''^k^{2gt,2gt')+Y. 

(D.12) 

with Km given by (14. 330 and 

g„(M-)-^°""'"];;!°"'"''". (D.13) 

Due to symmetry the second set of dressing roots can be treated in the same way. After 
elimination of the auxiliary densities in fID.Qp one finds 

2i9k{xi, . . . ,xm) = argcr^(a;i, . . . , xm) = 

.. / iti ., 

'29' 



/OO poo , / s 

dt e-^e**"'^ / dt' e-^ Yl [Kdi2gt', 2gt) - k,{2gt, 2gt')) (D.14) 
-oo J —oo i=l 

with the dressing kernel given by (II. 2p . Comparing the last equation with (ll.ip we notice 
that the dressing phase factorizes as promised in (ID.Tp . 

The above procedure is suggestive, but currently plagued by the following problems: 

• There is no solution to the equations (ID.2p . flD.5p relating to the dressing roots when 
their number Md is finite. Therefore expressions such as ( ID. lip , where we employ 
Md as an infrared regulator, are somewhat ill-defined. This is in contradistinction to 
our treatment of JF^ where the equations certainly make sense at finite L. 

• The scattering pattern of fID.ip - fID.Sp corresponds to a very curious "Dynkin di- 
agram" whose Lie-algebraic origin is very questionable, see Figure 2. Furthermore, 
the occupation (filling) numbers appear to be inconsistent with a standard nested 
Bethe ansatz. It would have to be shown how such filling numbers can be obtained 
from the non-compact nature of the underlying system. 

• The ansatz ( ID.ip - ( ID.SP leads to Bethe equations which are at Md = oo strictly 
identical to the original ones with the dressing factor. So the problem of the asymp- 
totic character of the equations of [Hj is not improved, and finite size effects are still 
not properly implemented. 
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® — o — ® — 6 — ® — o — ® 

Figure 2: Schematic scattering notation for the dressed asymptotic Bethe ansatz equations 
(1D.1I) - (ID.SIl . The horizontal branch represents the Dynkin diagram of psu(2,2|4) and 
the scattering of elementary magnons. The vertical branch (dotted lines) indicates the 
emulation of the dressing phase by the scattering of the dressing roots among themself 
and the momentum carrying main node. Note that the auxiliary "dressing" roots do not 
scatter with any of the auxiliary "nesting" roots. 

• Finally, it is completely unclear which "hidden" excitations are supposed to be scat- 
tering according to (]D.1|) - flD.51) . It is certainly not clear what, if anything, is being 
diagonalized by these equations. The equations might however turn out to be use- 
ful for proving the crossing invariance of the dressing factor at finite values of the 
coupling constant. 
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